Dynamics of a non homogeneously coarse grained system 
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To study materials phenomena simultaneously at various length scales, descriptions in which matter 
can be coarse grained to arbitrary levels, are necessary. Attempts to do this in the static regime (i.e. 
zero temperature) have already been developed. In this letter, we present an approach that leads 
to a dynamics for such coarse-grained models. This allows us to obtain temperature-dependent 
and transport properties. Renormalization group theory is used to create new local potentials 
model between nodes, within the approximation of local thermodynamical equilibrium. Assuming 
that these potentials give an averaged description of node dynamics, we calculate thermal and 
mechanical properties. If this method can be sufficiently generalized it may form the basis of a 
Molecular Dynamics method with time and spatial coarse- graining. 

PACS numbers: 83.10.-y, 62.20.-x, 81.40.Lm, 82.20.Wt, 63.10.+a 



Predicting macroscopic properties of materials start- 
ing from an atomistic or electronic level description can 
be a formidable task due to the many orders of magni- 
tude in length and time scale that need to be spanned. 
A characteristic of successful approaches to this problem 
is the systematic coarse-graining of less relevant degrees 
of freedom in order to obtain Hamiltonians that span 
larger length and time scales. For example, in first - 
principles thermodynamics of crystalline solids, which is 
one of the best developed examples of micro to macro 
bridging, electronic and vibrational excitations are inte- 
grated out in order to obtain lattice - model Hamilto- 
nians that describe the substitutional degree of freedom 
jij. Monte Carlo simulations can then be used to sim- 
ulate the kinetic evolution of the system or to obtain 
its thermodynamics state function. The reason first - 
principles thermodynamics has been so well developed is 
that it deals with extensive (averaged) quantities of ho- 
mogeneous materials and minor inhomogeneities in real 
materials, such as interfaces or dislocations, have a minor 
effect on the thermodynamic functions. Another extreme 
is the study of mechanical properties, such as plasticity, 
where the property of interest (i.e. plastic yield) is deter- 
mined by discrete events (slip of dislocations), but over a 
very large scale. This type of problems requires the use 
of inhomogeneous coarse - graining methods: atomistic 
- level resolution may be required near the key features 
in the material (e. g. dislocations or grain boundary) 
and lower resolution is needed in between, in order to 
make the problem computationally tractable. An im- 
portant step towards such a coupled multi - length scale 
description was taken with the development of the qua- 
sicontinuum method (QCM) in which the behaviour of 
groups of atoms (nodes) are treated with a finite ele- 
ment scheme j^] . The interaction between atoms is typ- 
ically calculated with empirical potentials. In its origi- 
nal form QCM is essentially a method that improves the 
boundary conditions on atomistic regions and allows the 
boundary conditions of different atomistic regions to in- 



teract. Since QCM consists of an optimization of the 
energy, no time or temperature phenomena are present. 
Introducing temperature into QCM can be done by us- 
ing potentials that incorporate the entropy due to lattice 
vibrations, for example in a local Einstein description 
P,0. This is conceptually similar to the coarse - grain- 
ing of vibrations in first - principles thermodynamics |^| 
and leads to structures that are free - energy minimized. 
The study of time/temperature dependent phenomena is 
more difficult and requires the development of a dynam- 
ics for a system with an inhomogeneous level of coarse 
graining. To our knowledge, no formal development of 
this problem exists for a more general class of interac- 
tions, though methods that couple Molecular dynamics 
to continuum description have been studied Q and Fi- 
nite Elements dynamics has been developed in the har- 
monic approximation without time rescaling [Q. Other 
coarse-graining techniques have been developed: effec- 
tive Langevin dynamics procedures have been proposed 
to describe relaxation of macroscopic degrees of freedom 
H , and transition state theories have been generalized to 
boost time evolution, by modifying the shape of the sur- 
face of contant energy (hyperdynamics) |)| or using the 
Onsager-Machlup action |n| . None of these methods has 
employed simultaneous space/time coarse-graining. 

In this letter we present a suggestion for the dynam- 
ical modeling of 2D systems that exist simultaneously 
at different level of coarse graining in time and space do- 
mains. Coarse graining requires both a scheme to remove 
atoms, and a prescription to define potentials between 
the remaining atoms. For the approach advocated in 
this letter atoms are integrated out through bond mov- 
ing, similar to the Migdal-Kadanoff approach in renor- 
malization group jl3|] . New potentials can be defined in 
various ways, but an important aspect is that the coarse- 
grained system ultimately evolves to the same equilib- 
rium state of the fully atomistic one. Hence, our cri- 
terium for defining new potentials when removing atoms 
is that the partial partition function of the system re- 



mains unchanged. We assess the validity of the model, 
by comparing the elastic, thermodynamic, and transport 
properties of a non-homogeneously coarse - grained and 
fully atomistic model. 

As a matter of introduction, the coarse-graining is first 
considered for a one-dimensional system. Let us consider 
a finite one - dimensional chain of atoms with mass m 
that interact through nearest neighbor pair potentials. 
Each particle has kinetic and potential energy and the 
total Hamiltonian is: H[{qi,pi}], where qi is the posi- 
tion coordinate and Pi is the momentum. It will be as- 
sumed that the gradient of temperature along the chain 
is small enough for local thermodynamic equilibrium to 
exist. One coarse graining step is defined as removing 
every second atom. The energy Htu(qi, qi+2) between 
nearest neighbors of the coarse system (2nd neighbors in 
the full atomic system) is defined so as to conserve the 
partition function (the subscripts (1) indicate the 1st step 
of renormalization): 

e -PH m to, 9i+a ) = 1 JJ dqi+ldpi+1 e -/J[Hte,« + 0+*(* +l ,« +a )] - 

Integrating the momenta leads to the the reduced for- 
mula: 

e -/3[V(i)fe,9i +2 ,T)+F Cl) (T,i+l)] = J dq . +ie -0[V(q i ,q i +i)+V(q i+ i,q i+ : 

where -Fm(T, i + 1) is an excess free energy that does 
not depend on the positions (ft, 5*4-2) < It contains the 
entropy of atom i + 1 lost in the renormalization step. 
It is necessary to keep track of this quantity to calculate 
properly the extensive thermodynamic quantities of the 
system. Finally, from the last equation, it is possible to 
extract an effective potential Vm(ft,ft4-2,T') (with har- 
monic potentials the linear chain model can be solved 
analytically). This potential is now temperature depen- 
dent 0. 

It is not obvious that this choice for the coarse-graining 
algorithm leads to the correct dynamics. However is is 
well know that the interaction so defined between parti- 
cles i and i + 2 is correct in the long time limit, i.e. when 
the motion of i + 1 is much faster that i and i + 2. Hence 
low frequency dynamics in the coarse system will likely 
be better represented than high frequencies. However, 
we feel that the ultimate justification for this approach 
should be evaluated on the basis of a comparison of prop- 
erties of the coarse and fully atomistic system. This is 
investigated in this letter. Hence this potential generates 
some dynamics of particles i and i + 2 by the averaged 
interaction of particle i + 1. In a first approximation, we 
consider F(i\(T, i + 1) to be independent from (ft,ft4-2), 
and the potential Vm to contain all the possible spatial 
dependencies of the remaining coordinates. 

As is typical in dynamic renormalization group theo- 
ries, integrating out degrees of freedom, leads to a time 



rescaling as tm = b z t where b = 2 is our scaling fac- 
tor and z is the dynamical exponent |T^| . This ex- 
ponent does not affect equilibrium properties and will 
determined later. Coarse graining has to conserve the 
total mass, so, we take, raf^ = m* + mj+i/2 and 
to|^ 2 = + ^i+i/2. To describe the potential we 

take a simple expansion around the minimum a (a = 3A 
atomistic lattice spacing and m, = 30 x neutron mass), 
V(ft, ft+i) = J2 n ~ 9i+i ~ a n a)' n . A symmetric po- 
tential has a n = 1, Vn. If we assume that the potential 
V(i) has the same functional shape of V , we get recursive 

relations Vry-i) Vu) (with Vm) = V) which creates 
higher order coarse-grained potentials: 

V {j )(qi,q i+b j,T) = J2k u)n (T)(q t - q i+bi - Va {j)n a) n . 

n 

Our model system interacts through a potential with 
second and fourth powers hi = 8.8 • 10 3 [K/A 2 ], — 
1.8 • 10 5 [K/A 4 ], which is either symmetric (pt2 = «4 = 
1), or non-symmetric («2 = 0.999, «4 = 1.0315), and 
we keep only the second and fourth power coefficients 

^(j)2' ^0)4' a (i)2' a U)i f° r ever y renormalized potential. 
To extend the method to a 2D lattice, we use the Migdal- 
Kadanoff moving bonds (fig-|l|) approximation |ll| , to re- 
move atoms. 
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FIG. 1. Triangular lattice with bond moving approxima- 
tion, and example of non-homogeneous coarse graining 

The model and assumptions are tested by compar- 
ing the results of molecular dynamics simulations on the 
original fully atomic system and the non-homogeneously 
coarse grained one. Regions with different coarse grain- 
ing are considered with different time-evolution, because 
of the dynamical exponent z in the time scaling. The 
method is tested on elastic properties, thermodynamic 
quantities, thermal expansion, and finally on its heat 
transport. We see these as successfully more stringent 
tests to pass. A 2D triangular lattice with 225 x 31 
atoms (6975 atoms) and symmetric potential, is stud- 
ied in response to static tensile/shear stress, and under 
isotropic pressure deformation. The non-homogeneously 
coarse grained system is represented by 1510 nodes in 
an arrangement similar to that in figure [l[ We calcu- 
late the strain response to a normalized tensile stress 
{<j* x = <7\i/o-q) along the longest direction (rjo is the 
tensile stress which gives 5% of lattice distortion) . Vari- 



ous simulations are done at different normalized temper- 
atures ( T* = T/Tq = 0.04- • -0.64, where T = huj Q /k B , 
cjo = \fk-2jm). Figure || shows the strain response for 
T* = 0.32. The linear / non-linear regions of the strains, 
are clearly conserved by the non-homogeneous coarse 
graining. The elastic modulus of the two systems are 
equal to within 4%. Similar results have been obtained 
for shear strain (5%). For the bulk modulus calcula- 
tions we use the non-symmetric potential, as described 
in the previous section. In the range of temperatures 
considered the bulk modulus Bh for the original system 
(subscript h) and B n h for the non- homogeneous coarse 
grained one (subscript n h) are equal to within ±1% The 
good agreement for the elastic properties may not be sur- 
prising but indicates that the bond- folding does not mod- 
ify the macroscopic energies. As the elastic properties are 
largely a reflection of the direct interaction and are not 
much influenced by temperatures or atomic motion, they 
do not really test the assumption made on the coarse- 
grained dynamics. 
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a normalized stress 

FIG. 2. Strain versus normalized stress <x*i. Contin- 
uous line represents the strain e/ in for the original sys- 
tem, while squares □ represent the strain e n hu for the 
non-homogeneously coarse grained one (T* = 0.32) 

We also evaluate the heat capacity for the 2D system 
with non-symmetric potential. The heat capacity from 
the coarse system can not be directly compared to the 
complete system, but needs to be augmented with the 
contribution from the entropy that is lost by removing 
degrees of freedom. This entropy can be calculated by 
keeping track of all the free energies Fu) (T, i) produced 
by the renormalization integration and taking tempera- 
ture derivatives to get entropies Q] . We find that within 
the range of simulations T* = 0.1 • • • 1.6, Cy nh and Cvh 
are equal to within the numerical noise of ±0.2%. Since 
the calculation is classical, the specific heats do not go to 
zero for T* -> 0. 

The calculation of properties at constant pressure, 
such as thermal expansion or C p , requires particular ef- 
fort to reproduce the correct dependence of free energy 
on temperature and volume. This dependence is hid- 
den in either the nonlinear terms of the renormalized 
potentials or the volume dependence of the integrated 
free energy F(j\(T,i). An approximate representation of 
the thermal expansion can be achieved by interpolation 



of ay) ,Q!(j) to conserve the thermalized bond length 

< Qi+bJ - ft >(i)= b < 1t+{b-iy ~ ft >(j'-i)' ( the sub " 
scripts on the averages indicate the potential used for 
the calculation). This constrained can be satisfied with 
a very small shift (« 10~ 5 ~ 10~ 4 ) of the parameters 
<2(j) 2 ! a (j)4- The overall effect is simply an effective bias 
on the potential to match the correct thermal expansion. 
Simulations with this method for the thermal expan- 
sion coincide for the homogeneous and the non homoge- 
neous coarse grained systems, within the numerical noise 
(±5%). It is important to understand that without this 
correction, coarse-grained systems have incorrect thermal 
expansion (due to the removal of volume-dependent en- 
tropy). Hence simulations with inhomogeneously coarse 
grained regions would build up large internal strains upon 
changing the temperature, a fact that does not seem to 
have been recognized in previous formulations. 

As stated in the description of the model, it is nec- 
essary to calculate a dynamical property (such as the 
thermal conductivity k), to get the value of the dy- 
namical exponent z. For finite systems, we can com- 
pare the results for a homogeneous lattice to those of a 
non-homogeneously renormalized one, if we assume the 
effects of finite size to be the same. To determine k, 
several simulations for a 2D homogeneous triangular lat- 
tice (N x x N y — 225 x 33) and for its homogeneously 
coarse grained equivalent are run at various tempera- 
tures. A non-symmetric potential which allows anhar- 
monic phonon interactions, is used. Since there are no in- 
terfaces in the homogeneously coarse grained lattice, the 
only effect is due to the time scaling. The thermal con- 
ductivity k is obtained from the standard Green-Kubo 
relation p5[ . This relation involves time integration, so 
the dynamical exponent z can be implicitly obtained by 
fitting k which is defined per "unit time". We find that 
for a value of z w 1.45 ± 0.1, perfect agreement exists for 
the k in the atomistic and coarse system. 
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FIG. 3. Interface effect: thermal conductivity ratio 
K nh (N™ h )/K h (N^) versus the length of a 2D lattice with one 
interface. (+) for T* = 1 and (o) for T* = 3. 

In order to study the thermal transport in non- 
homogeneous systems the effect of interfaces between dif- 
ferently coarse-grained regions need to be understood. 
We run simulations for a 2D homogeneous lattice, with 
N h = N£ x atoms, and for its non homogeneous 



coarse grained version, with N nh — N™ h x iV™' 1 nodes. 
The non homogeneous lattice has one coarse-graining in- 
terface at the center parallel to y direction. To one side 
of this interface everything is fully atomistic, while to 
the other side one level of coarse-graining is applied. 
Thus N£ h = 2N%/3, and each region is A^/3 wide. 
By systematically varying the size of the systems and 
plotting K nh (N£ h )/n h (N£), we determine the effect of 
the interface/size in the thermal conductivity. Figure 
H shows this effect at T* = 1 and 3: at high N%, the 
length of the regions becomes comparable to the phonon 
mean free path. Hence the fact that there are less nodes 
in the coarse-grained system becomes less apparent and 
n nh (N" h ) -> n h {N^). The picture shows that, for our 
potential, dynamical exponent z = 1.45, and temper- 
atures (T* = 1 and 3), regions wider than 500 nodes 
give acceptable results. 
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FIG. 4. Fhermal conductivities K h (T*)/ k h (T* = 1) (o) and 
K ni>(T*)/ K h (T* = l) (+) of 2D lattices for various normalized 
temperatures. 

For high enough separation between interfaces, the 
coarse-grained dynamics reproduces the T-dependence 
of the thermal conductivity well. Figure ^ shows 
K nh < h (T*)/n h (T* = 1) for a 2D lattice with N h = 
N% x N' y l = 11264 x 33 = 371712 atoms, and its non- 
homogeneous coarse grained version with N — 63072 
nodes, and 6 interfaces parallel to y axis (each region 
is 512 nodes wide at different level of coarse graining). 
The non-homogeneous model underestimates the ther- 
mal conductivity at low temperatures (T* < 1/3). At 
high temperatures (T* > 1) the two results are within 
10 ~ 15%. Normal materials have To of the order of 
room temperature, so the error is acceptable considering 
the strong approximations in the method. 

To conclude, we have proposed and analyzed a molec- 
ular dynamics method with time and spatial coarse- 
graining. We show that mechanical and thermodynam- 
ical properties are in excellent agreement with the non 
coarse grained system. If this method can be sufficiently 
generalized in 3D, it may form the basis of a RG Molec- 
ular Dynamics to investigate effects of temperature and 
defects in real nanostructures. 
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